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Abstract 

In this paper Quintic Spline is defined for the numerical solutions of the 
fourth order linear special case Boundary Value Problems. End condi- 
tions are also derived to complete the definition of spline. The algorithm 
developed approximates the solutions, and their higher order deriva- 
tives of differential equations. Numerical illustrations are tabulated to 
demonstrate the practical usefulness of method. 



1 Introduction 

Spline functions are used in many areas such as interpolation, data fitting, numerical solu- 
tion of ordinary and partial differential equations. Spline functions are also used in curve 
and surface designing. 

Riaz A. Usmani [1], considered the fourth order boundary value problem to be the problem 
of bending a rectangular clamped beam of length 1 resting on an elastic foundation. The 
vertical deflection w of the beam satisfies the system 

d 4 

w = D 1 q(x) where L = -— (1.1) 

ax 4 

w{0) = w(L) =w'(0) =w'(L) = (1-2) 

where D is the flexural rigidity of the beam, and k is the spring constant of the elastic 
foundation and the load q(x) acts vertically downwards per unit length of the beam. The 
detail of the mechanical interpretation of (1.1) belongs to a general class of boundary value 
problems of the form 




2 



d 



•A 



dx 4 + f(x))y(x) = g(x), x G [a, b] 
y(a) = a , y{b) = a x 
y'(a) = /?o, y'(b) = /9i 



(1. 3 



where ctj, (3i ; i= 0,1 are finite real constants and the functions f(x) and g(x) are continuous 
on [a,b]. The analytic solution of (1.3) for special choices of f(x) and g(x) are easily obtained, 
but for arbitrary choices, the analytic solution cannot be determined. 
Numerical methods for obtaining an approximation to y(x) are introduced. Usmani [1] 
derived numerical techniques of order 2, 4 and 6 for solution of a fourth order linear boundary 
value problem. Usmani [2] derived cubic, quartic, quintic and sextic spline solution of 
nonlinear boundary value problems. Usmani and Manabu Sakai [3] developed a quartic 
spline for the approximation of the solution of certain two point boundary value problems 
involving third order linear differential equation. 

N.Papamichael and A.J.Worsey [4] derived end conditions for cubic spline interpolation 
at equally spaced knots. N.Papamichael and A.J.Worsey [5] have developed a cubic spline 
method, similar to that proposed by Daniel and Swartz [6] for second order problems. In this 
paper a quintic spline method is described for the solution of (1.3). The end conditions for 
quintic spline interpolation, at equally spaced knots are derived which uniformly converge 
on [a,b] to 0(/i 6 ), which is discussed in the next section. 



2 Quintic Spline 

Let Q be a quintic spline defined on [a, b] with equally spaced knots 

Xi = a + ih; i = 0, 1,2, k ( 2. 1 ) 

where 

* = ^i (2.2) 

Moreover for i=0,l,2,. . . , k, taking 

Q{ Xi )=y t , Q^( Xi ) = m h ( 2. 3 ) 

Q^(xi) =M t , Q®(xi) = m, ( 2. 4 ) 

and 

Q^( Xi ) = Ni (2.5) 

Also, let y(x) be the exact solution of the system (1.3) and y\ be an approximation to y(xi), 
obtained by the quintic spline Q{xi). It may be noted that the Qi(x), i = 1, 2, 3, . . . , k can 
be defined on the interval [xj_i,Xj], integrating 

Q\{x) = ^[JVi_i(xi - x) + N t (x - xi_i)] ( 2. 6 ) 

four times w.r.t. x, which gives 

1 /It 3 /?t 2 

Qi{x) = —[N i _ 1 {x i -xf + Niix-Xi-i) 5 ] + — + — + Cx + D (2. 7) 
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To calculate the constants of integrations, the following conditions are used. 

Qi{xi) = yi ; Q 2 i(xi) = Mi 
Qi{xi-i) = yi-i ; Q 2 i(xi„i) = ( 2. 8 ) 

The identities of quintic splines for the solution of (1.3) can be written as 

5 . 

mi-2 + 26mj_i + 66mj + 26m i+ i + m i+2 = - [-yi-2 - 10yi-i + 10yi+i + yi+2\ , 

* = 2,3, As — 2 ( 2. 9 ) 

20 

Mi_ 2 + 26Mi_! + 66Mi + 26M i+1 + M i+2 = + 2y*-i - % + 2y m + y i+2 ] , 

i = 2,3,...,fc - 2 ( 2. 10 ) 

60 

rii_ 2 + 26nj_i + 66rii + 26n m + n i+2 = [-yi-2 + - 2y i+ i + yi+ 2 ] , 

i = 2,3,...,fc-2 ( 2. 11 ) 

120 

iVi_ 2 + 26iVi_i + 66iVi + 26N i+1 + N i+2 = — [ yi _ 2 - 4y;_i + 6^ - % +1 + y i+2 ] , 

i = 2,3,...,fc-2 (2. 12 ) 

JVi_i + ANi + iV i+1 - p[A4i_i - 2Mi + M i+ i] = , 

i = 1,2,3,..., k - 1 ( 2. 13 ) 

60/i{mj_i + 2mj + m i+1 } - /i 3 {3raj_i + 14nj + 3ra m } = 120{-yj_i + y i+ i} , 

i = 1,2,3,. ..,k - 1 ( 2. 14 ) 

8/i{m m - - /i 2 {Mi_i - 6M; + M m } = 20{y^i - 2 Vi + y i+1 } , 

i = 1,2,3,..., k - 1 ( 2. 15 ) 

h h 3 1 

mi = -{2Mi + Mi.x} -—{8Ni + 7iVi_i} + -{y;-y;_i} , 

i = 1,2,3, ...,k ( 2. 16 ) 

= -J{2Mi + M^ + ^-SN, + 7iV m + \{y i+ i - Vi } , 
o 360 h 

i = 0,1,,..., k - 1 ( 2. 17 ) 



i = 1,2,3,...,*; - 1 ( 2. 18 ) 
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h 2 1 

Mi = ~J^{Ni-i + 8N t + N i+1 } + j^{yi-i ~ 2Vi + Vi+i} , 

i = 1,2,3, ...,k - 1 ( 2. 19 



Mi = 7^{ m i-2 + 32mi_i - 32m m - m i+2 } 
5 r 

+ {2/i-2 + 16yi_i - 34y; + 16y i+ i + y; +2 } , 

i = 2,3,..., A; - 2 ( 2. 20 ) 

= {Mi_!+ 18Mi + M m } + $ {2/i_i - 2 Vi + y m } , 

i = 1, 2, 3, k - 1 ( 2. 21 ) 



The relations (2.9)— (2.21) can be derived from the results of Albasiny and Hoskins(1971), 
Fyfe(1971), Ahlberg, Nilson and Walsh(1966) and Sakai(1970) discussed by N. Papamichael 
[7]. The uniqueness of Q can be established showing that any of the four (k + 1) x (k + 1) 
linear systems, obtained, using one of the relations (2.9), (2.10), (2.11) and (2.12)together 
with the four end conditions. In this paper the linear system corresponding to (2.12) has 
been chosen and has a unique solution for NiS , i = 0, 1, k. Equation (2.13) and (2.19), 
or (2.19) and (2.21) give the parameters MjS , i = 0,1, k. Consider the system (2.12) 

120 

N^ 2 + 26iVi_i + 66iVi + 26iV m + N i+2 = [y^ 2 - 4^_i + 6 Vi - 4y i+1 + y i+2 } , 

i = 2,3,...,k -2 ( 2. 22 ) 

where 

Ni = -fiyi + g t ( 2. 23 ) 

The above system gives (k — 3) linear algebraic equations in the (k — 1) unknowns (yi, i = 
1, 2, k — 1). Two more equations are needed to have complete solution of yis which are 
derived in section 3. 

It may be noted that Papamichael [5] needed two consistency systems for the 
solution of boundary value problem (1.3) but the method developed in this 
paper, needs only one such system. 

Taking forward difference operator E = e hD , equation (2.12) can be rewritten as 

120 



[E~ 2 + 2§E~ X + 66/ + 2§E X + E 2 )Ni = —j- [ E~ 2 - 4E" 1 + 61 - AE 1 + E 2 } y t , 

i = 2,3,..., k -2 ( 2. 24 ) 



AT. = ^ ( e - 2hD + 26 e- hD + 66 I + 26 e hD + e 2hD y l 
{e -2hD _ 4 e -hD + g j _ 4 e hD + e 2h D) y . ) 

i = 2,3, ...,k — 2 ( 2. 25 ) 
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Expanding the R.H.S. of above, in terms of power series and dividing give, 



12 



240 



( 2. 26 ) 



Papamichael [7] proved the following lemma to determine the end conditions in terms of 
first derivative of quintic spline 



Lemma 1 

Let Aj = rrii — If y G C 7 [a, 6] then 



Ai_ 2 + 26Ai_i + 66Ai + 26A m + X i+2 = ft , i = 2, 3, .., k — 2 



( 2. 27 ) 



where 



< ^ /i 6 ||y (7) || , i = 2,3,...,fc-2 

Using the above lemma along with equation(2.12) and with Taylor series expansion about 
the point Xj, the following lemma can easily be proved, for the case discussed in this paper. 



Lemma 2 

Let 



Al ~ Vi 12 Vi + 240 y * 



If y € C 10 [a,6] then 

Ai_ 2 + 26Ai_i + 66Ai + 26A m + \ i+2 = ft , i = 2, 3, fc - 2 

where 

2665920 



( 2. 28 ) 
( 2. 29 ) 



IAI < 



10! 



/i 6 ||y( 10 )|| , i = 2,3,...,fc-2 



Finally the required end conditions are derived in the following section. 



3 End Conditions 

Consider end conditions of the form 

A?o + aN x + 0N 2 + <yN 3 + A/4 
and 



/i 4 



^aiyi + bhy^ + h A ^qy, 



.(4) 



L i=0 



i=0 



(3. 1) 



N k + aN k ^ + /3Af fc _ 2 + 7 iV fc _ 3 + N k _ 4 = 



3 4 
i=0 i=0 

( 3. 2 
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where the scalars a, (3, 7 , cii, i=0,l,2,3, b and q, i=0,l,...,4 can be determined with the 
assumption that Q exists uniquely and 



\\Q r - 2/1 = 0(h e ~ r ) , r = 0,l,...,5 

For this let 

\ . _ ? / 4 ) h * „(6) , ^ ..(8) at. ,-_ n 1 fr 

Since it is supposed that y € C 10 [a, 6], therefore the equations (2.12), (3.1) and (3.2) give, 

A + aAi + /3A 2 + 7 A3 + A 4 = Pi ( 3. 3 ) 



Ai-2 + 



where 
Pi = 



26Ai_i + 66Ai + 26A m + X i+2 = fa , 

i = 2,3, ...,k - 2 ( 3. 4 ) 

Afe + aA fc _i + /?A fc _ 2 + 7^-3 + A fc _ 4 = /3 fc _i ( 3. 5 ) 



. i=0 i=0 



and 

Pk-i 



1 

/i 6 , (6) (( 



■ (4) + y?> 



(3. 6) 



1 



-J2am- t - bhyM - + (y W + ^ + ^(4) 2 

i=n 



i=0 

(4) (4) x h 6 ( 6 ) ( 6 ) ( 6 ) ( 6 ) ( 6 ) 

+ 72/fc_ 3 + 2/fc- 4 ) " ^ ( i/jfe + «2/fc-l + + 71/fc-3 + 24-4 

„(8) 



and, from lemma 2 



2665920 , fi „ nm „ 
l^l < 1Q , h 6 \\y^\\ , * = 2,3,..,fc-2 

Following Papamichael [7], the required end conditions may be written as 
.... 1 r 220 



(3. 7) 



(3. 8) 



iV + Ni = ^_ _ I20yi + 60y 2 - yy 3 + 40/iyJ ; 

4 f 2519 ( 4 ) , 11800 00Q 



+ /,4 f 2519 „(4) , 11822 „(4) , 22 3 (4) 176 

+ ^ V 2520 2/0 + "l260~ 2/1 + 210 2/2 + 180 



1769 (4) ^ 
2520 tf4 , 



(3. 9) 
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1 



220 40 

— y k - 120y fc _i + 60y fc _ 2 - — y fe _ 3 + AOhy^ 



and 

N k + iV fc _i 



( 3. 10 

The quintic spline solution of the system (1.3) is defined in the next section. 



+ n 1 2520 y * + 1260 Vk ~ l + 210 + 180 Vk ~ 3 

1769 ( 4 ) 



+ 2520 Vk ~ A 



4 Quintic Spline Solution 

The quintic spline solution of (1.3) is based on the linear equations (2.12), (3.9) and (3.10). 
Let Y = yi, C = Cj, e = e*. Then the parameters j/j of Q satisfy the following matrix 
equation 

(A + h 4 BF)Y = C + e 

where Y, C, e are (k — 1) dimensional column vectors and A, B, F are (k — 1) x (k — 1) 
matrices, where 
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r 31686 


669 


528 


5307 




5040 


8400 


7200 


100800 




13 


11 


13 


1 




60 


20 


60 


120 




1 


13 


11 


13 


1 


120 


60 


20 


60 


120 



B 



1 


13 


11 


13 


1 


120 


60 


20 


60 


120 




1 


13 


11 


13 




120 


60 


20 


60 




5307 


528 


669 


31686 




100800 


7200 


8400 


5040 



F = diag(fi) 



11 OI (l) /i 4 /-3 3 , 31686 



669 



3696 5307 

53 + -T^TT 94 



180 



360 



° 2 = 120 90 ~ ^° V ° + 26 91 + 66 52 + 26 53 + 54 



yo 



(4. 1) 
(4. 2) 



Y^Tj ( 3i-2 + 26 gi-i + 66 & + 26 g i+1 + ) 
i = 3, 4, fc — 3 



(4. 3) 



Cfc-2 



120 



( 9k ~ fk Vk + 26 + 66 # fc _ 2 + 26 g k _ 3 + c/ fe _ 4 ) - y k ( 4. 4 ) 



(i) /i 4 /-3 3 

I + Qk + 

280 V360 360 



31686 

9k + tt^t: /fc ^ — tt^t ffe-i 



180 



669 3696 5307 

+ -gQ- fffc-2 + 18Q fffc-3 + 360 



(4. 5) 



and ej = 0(/i 6 ) i = 1, 1 

Extending the method for the solution of sixth, eighth and higher order bound- 
ary value problems, is in process. 

To implement the method for the quintic spline solution of the boundary value problem 
(1.3), three examples are discussed in the following section. 
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5 Numerical Examples 

In this section numerical technique discussed in section 4 is illustrated, by the following 
three boundary value problems of the type (1.3). 

Example 1 

Consider the following boundary value problem 

y iv + Ay = 1, y(-l) = y(l) = 0, 1 

( 5. 1 ) 

i ( n _ _ sinh(2) - sin(2) i 

y y~ L > - - y y L > - 4(cosh(2) + cos(2)) J 

The analytic solution of the above problem is 

y{x) = 0.25[1 — 2[sin(l) sinh(l) sin(z) sinh(x) 

+ cos(l) cosh(l) cos{x) cosh(x)]/(cos(2) + cosh(2))] 

The results are summarized in Table 1. 

Table 1: Maximum absolute errors for Problem 5.1 in y\^\ /i = 0, 1,2,3,4 



h 


yi 




(2) 

y\ 


(3) 


(4) 

Vi 


1 

8 


2.1 x 1(T 3 


5.5 x 10~ 3 


9.7 x 10~ 3 


1.06 x 10~ 2 


8.5 x 10~ 3 


1 

16 


9.10 x 1(T 5 


3.34 x 10" 4 


8.86 x 10" 4 


1.4 x 10~ 3 


3.64 x 10~ 4 


1 

32 


5.12 x 1(T 6 


1.77 x 10~ 5 


1.20 x 10~ 4 


2.60 x 10~ 4 


2.05 x 10~ 5 


1 

64 


2.85 x 1(T 6 


4.50 x 10~ 6 


2.44 x 10' 5 


6.03 x 10~ 5 


1.14 x 10~ 5 


1 

128 


7.87 x 1(T 7 


1.18 x 10" 6 


5.83 x 10~ 6 


1.47 x 10~ 5 


3.23 x 10~ 6 


1 

256 


1.98 x 1(T 7 


2.79 x 10~ 7 


1.44 x 10~ 6 


3.59 x 10~ 6 


8.73 x 10~ 7 


1 

512 


4.15 x 1(T 8 


4.04 x 10~ 8 


3.05 x 10~ 7 


6.81 x 10~ 7 


2.46 x 10~ 7 


1 

1024 


1.07 x 10~ 7 


1.91 x 10~ 7 


8.01 x 10~ 7 


2.12 x 10~ 6 


3.49 x 10~ 7 
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Example 2 

Consider the following boundary value problem 

( 5. 2 



y ™ + xy = _(8 + 7x + x^e 21 , 
y(0) = y(l) = 0, y'(0) = 1, y'(l) = -e 



The analytic solution of the above differential system is 

y(x) = x(l — x) e x 

The observed maximum errors (in absolute value) associated with yf^, n = 0, 1, 2, 3,4 , for 
the system 5.2 are briefly summarized in Table 2. 



Table 2: Maximum absolute errors for Problem 5.2 in y\ , /i = 0, 1, 2, 3, 4 



h 


Vi 




(2) 

Vi 


(3) 

Vi 


(4) 

Vi 


1 

8 


5.3828 x 10~ 4 


1.9 x 10~ 3 


6.0 x 10~ 3 


3.8 x 10~ 2 


3.2768 x 10~ 4 


1 

16 


8.8114 x 10~ 5 


2.9045 x 10~ 4 


1.4 x 10~ 3 


9.9 x 10~ 3 


5.3779 x 10~ 5 


1 

32 


1.636 x 10~ 5 


5.1554 x 10~ 5 


3.3586 x 10- 4 


2.6 x 10~ 3 


9.824 x 10~ 6 


1 

64 


3.3894 x 10~ 6 


1.058 x 10~ 5 


8.2872 x 10~ 5 


6.6254 x 10~ 4 


1.991 x 10~ 6 


1 
128 


7.5919 x 10~ 7 


2.388 x 10~ 6 


2.269 x 10" 5 


1.6824 x 10~ 4 


4.3958 x 10~ 7 


1 
256 


1.7884 x 10~ 7 


5.671 x 10~ 7 


5.9302 x 10~ 6 


4.2401 x 10~ 5 


1.0264 x 10~ 7 


1 
512 


3.7364 x 10~ 8 


1.1843 x 10~ 7 


1.3068 x 10~ 6 


9.4838 x 10~ 6 


2.1304 x 10~ 8 


1 

1024 


6.3664 x 10~ 8 


2.001 x 10~ 7 


1.9244 x 10~ 6 


9.969 x 10" 6 


3.5439 x 10~ 8 



11 



Example 3 

Consider the differential system 

U lv -y = - 4(2x cos (x) +3sin(x)), 1 , , 

1/(0) = = , y'(-l) = 2 sin(l) , y'(l) = 2 sin(l) J ( 5 ' j 

The analytic solution of the above system is 

y(x) = (x 2 — 1) sin(x) 

The observed maximum errors (in absolute value) associated with yf\ n = 0, 1, 2, 3, 4 , for 
the system 5.3 are briefly summarized in Table 3. 



Table 3: Maximum absolute errors for Problem 5.3 in y\ , fi = 0, 1, 2, 3, 4 



h 


Vi 




(2) 

Vi 


(3) 

Vi 


(4) 

Vi 


1 

8 


1.4 x 1(T 3 


9.9 x 10~ 3 


4.3 x 10~ 2 


1.0 x 10" 1 


1.4 x 10" 3 


1 

16 


1.42 x 10~ 4 


1.2 x 10~ 3 


7.0 x 10 -3 


2.26 x 10~ 2 


1.42 x 10~ 4 


1 

32 


8.18 x 1(T 6 


1.34 x 10~ 7 


1.2 x 10~ 3 


5.1 x 10~ 3 


8.81 x 10~ 6 


1 

64 


5.05 x 10~ 6 


1.95 x 10~ 5 


2.51 x 10~ 4 


1.2 x 10~ 3 


5.05 x 10~ 6 


1 

128 


1.6 x 10~ 6 


5.88 x 10~ 6 


5.68 x 10~ 5 


3.01 x 10" 4 


1.61 x 10~ 6 


1 

256 


4.42 x 10~ 7 


1.58 x 10" 6 


1.35 x 10^ 5 


7.48 x 10~ 5 


4.42 x 10~ 7 


1 

512 


1.15 x 10~ 7 


4.09 x 10~ 7 


3.31 x 10~ 6 


1.86 x 10~ 5 


1.15 x 10~ 7 


1 

1 024 


3.19 x 10~ 8 


1.05 x 10~ 7 


8.51 x 10~ 7 


5.01 x 10~ 6 


3.19 x 10~ 8 
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